Anomalous temperature dependence of the the spin-flip thermalization time between 
the dark and bright exciton states in silicon nanocrystals 
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Silicon nanocrystals are studied by time-resolved fluorescence spectroscopy. After laser excitation 
the bright and dark exciton ground state levels are populated at random, but subsequently the 
decay curves reveal a thermalization between these levels. The characteristic thermalization time 
is found to be approximately 100 ns for temperatures below 100 K and surprisingly increases for 
higher temperatures. The decay curves are analyzed using a simple two- state model for the bright 
and dark exciton ground states. 



Since the discovery of light emission from porous 
silicon 1 , the structural and optical properties of nano- 
structured silicon has been studied extensively. In par- 
ticular, the electron- hole exchange interaction^ is respon- 
sible for splitting the exciton ground state into bright (ra- 
diative recombination dipole allowed) and dark (almost 
dipole forbidden) states, which was demonstrated-^^ 
by time-resolved fluorescence measurements for various 
temperatures while assuming the bright and dark state 
populations to be in thermal equilibrium. However, the 
spin-flip mechanism behind this thermalization has not 
yet been studied for silicon nanocrystals (NCs). In this 
letter we demonstrate that the spin-flip process can be 
seen directly in the luminescence decay curves, and we 
measure the characteristic time scale for the thermaliza- 
tion versus temperature. 

A sample of silicon NCs was prepared by magnetron 
sputtering, annealed at 1100°C for 1 hour in N2 (2 Bar), 
and subsequently passivated at 500° C for 1 hour in 95% 
N 2 + 5% H 2 (2.4 Bar). By co-sputtering Si and Si0 2 on a 
Si- wafer an approximately 300 nm thick layer of SiO^ was 
achieved with x = 1.601q! !io determined by Rutherford 
backscattering spectrometry (corresponding to a silicon 
excess concentration of 25^%). The room-temperature 
luminescence spectrum (late-time spectra in Fig. [2(c)) 
is peaked at 787 nm, being consistent with previously 
reported results for ~ 2 nm diameter NCs 7 (when con- 
sidering the differences^ between magnetron sputtering 
and plasma-enhanced chemical vapor deposition). 

The experimental setup is shown in Fig. [2(a). A 
frequency-doubled Ti:sapphire femtosecond laser deliv- 
ers pulses at 400 nm with a repetition rate of 1 kHz and 
excites the NC sample, which is contained in a closed- 
cycle cryostat. The fluorescence is spectrally filtered by 
a monochromator (AA = 2.5 nm) and detected by a sil- 
icon avalanche photo-diode (APD). A reference detector 
enables correction for variations in pump power. The 
sample emits light in two main bands, see Fig. [2(c). At 
early times the luminescence is centered around 550 nm 
and decays on a time scale of 1.5 ns independent on emis- 
sion wavelength and temperature. For this reason, the 
long-wavelength tail of the early-time band can be iden- 
tified and subtracted from the late-time band, which is 



centered in the range of 750-800 nm. We attribute the 
late-time emission to the exciton recombination in the 
NCs 3 while the early-time luminescence closely resem- 
bles the emission from oxygen-related defects 9 . Support 
to this picture is also given by the fact that the late-time 
emission center energy varies with temperature in a way 
not very different from the bulk silicon bandgap variation 
(Fig. [2(d)), while the early-time band center energy is in- 
dependent on temperature. In the following we focus on 
the late-time band only, and the early-time contribution 
has been subtracted from all the data shown in Fig. [2j 
The ground state of the quantum confined exciton is 
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FIG. 1: (a) Experimental setup, (b) The exciton model with 
spin- flip and decay rates defined. |B) and |D) denote the 
bright and dark exciton state, respectively, while |g) is the 
crystal ground state, "rad" and "nrad" denote radiative and 
non-radiative decay channels, (c) Time-gated spectra. The 
early-time spectra are detected for t < 25 ns and shown in 
red (294 K) and magenta (16 K). The late-time spectra are 
detected for t > 100 ns and shown in blue (16 K) and green 
(294 K). (d) The shift in peak center position of the early- 
time (blue) and late-time (red) spectra versus temperature. 
The bandgap energy shift of bulk is shown in black. 



split into two (double-degenerate) states, a bright state, 
|B), and a dark state, |D) (Fig. [TJb)). The bright 
state can recombine radiatively via a A J = 1 transi- 
tion, while the dark state is ideally radiatively forbidden 
(A J = 2) 2,4 . The selection rules do not apply strictly, 
but we expect a small ratio, R = rD,rad/rB,rad, between 
the radiative decay channels from the exciton states. 

The simple model of Fig. [2(b) has previously 
been applied successfully to other semiconductor 
NC a 10 ! 11 ! 12 ! 13 ! 14 ! 15 and is described by the rate equations: 
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where pB and describe the populations in the bright 
and dark exciton states, respectively. The spin-flip rates 
between |B) and |D) are denoted by Tbd and Tdb- 
Tb = TB,rad + rB,nrad is the total decay rate of the bright 
state (and similarly for Td). Equations (pQ) can be solved 
exactly, leading to a bi-exponential decay of pB and pB>. 
However, before writing the solution, we make the fol- 
lowing simplifying assumptions: (1) The exciton spin- flip 
rates are fast compared to the exciton decay rates such 
thatpife T B d, r DB > r B ,r D . (2) The two rates T B d and 
Tdb must drive the populations pB and toward ther- 
mal equilibrium. The relation Tdb = TBD^~ A ^ kT will 
assure this, where A is the energy splitting between the 
bright and dark states, k is Boltzmann's constant, and 
T is the temperature. (3) Immediately after laser excita- 
tion and subsequent carrier relaxation (which takes place 
on a sub-picosecond time scaled) to the exciton ground 
states the population of the exciton states is random: 
Pb(0) = p D (0) = \. We then get: 
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where the fast decay rate characteristic of the spin ther- 
malization is given by: 7f = Tbd + Tdb, and the slow 
decay rate, 7 S = rB 1 e +e _ A/A ^ rD , characterizes the total 

population decay: Pb(^) + pB>(f) = e _7s£ . Note that for 
times, t ^> 7 f _1 , the exciton populations have thermal- 
ized: p B (t)/p B (t)=e- A / kT . 

The time-dependent probability of photon emission de- 
pends on the radiative decay rates: p(t) = TB,r&dpB(t) + 
TD,radPD(£), and is written: 
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We stress that at zero time, p(0) = - 9 
pendent on temperature. This reflects the initial ran- 
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FIG. 2: Decay curves for temperatures 294 K (red), 100 K 
(green), 60 K (blue), and 16 K (black). The vertical axis 
is the same for all curves in (a) and (b). (a) The initial 
fast thermalization of fluorescence, (b) The long-time decay 
curves. The red dashed line denotes the level of dark counts, 
(c) The 16 K decay curve normalized to the black line of panel 
(a). The blue line is a bi-exponential fit. (d) As panel (c) but 
with the 294 K data. The blue line is a single-exponential fit. 



dom population of the exciton states. At longer times, 
t ^> 7 f _1 , the second term determines the light emission. 
Experimentally, the detected fluorescence from a sample 
containing NCs will depend on NC density, excitation 
power, detection efficiency of the entire optical setup, 
etc. Hence an unknown front factor must be multiplied to 
Eq. (j4j) and we cannot determine TB,rad and rD,rad on an 
absolute scale but only the ratio, R = r D ,rad/rB ; rad? as 
was also pointed out previously-^^. However, the ratio, 
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, is independent on the specific 



Af 

A s ~ R+e -A/kT 

experimental setup. Information about the total decay 
rates, Tb and Td, and the spin- flip rates, Tbd, Tdb, can 
be extracted from 7 S and 7f, respectively. In practice, 
the decay curves will not be single-exponential due to in- 
homogeneous broadening of the decay rates. However, it 
can easily be shown that the predictions for Af and A s 
are valid provided that TB, r ad and rD,rad represent the 
mean values of the radiative decay rates. 

Time-resolved decay curves were obtained for nine dif- 
ferent temperatures between 16 K and 294 K at detection 
wavelengths following the center of the late-time emis- 
sion spectra (Fig. [2(c)). Four representative curves are 
shown in Fig. [2](a,b). In Fig.[2](a) it can be seen that the 
zero-time fluorescence is essentially independent on tem- 
perature while it is very different after one microsecond. 
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FIG. 3: (a) The measured value of A s versus temperature. 
The vertical axis is normalized to the 16 K data point, (b) 
The measured ratio, Af/A a , versus temperature. In panels 
(a) and (b) the red curve follows the model of Eq. (U} with 
A = 15.0 ± 1.5 meV and R = 0.013 ± 0.003. (c) The charac- 
teristic decay time of the exciton levels versus temperature, 
(d) The characteristic time of the bright/dark-exciton state 
thermalization versus temperature. 



This is consistent with our assumption of initial random 
population in |B) and |D) followed by thermalization. In 
Fig. [2](b) the decay curves are shown for the entire laser 
repetition period of 1 ms. The characteristic decay time 
decreases with decreasing temperature since the popula- 
tion freezes out£^££ in the dark state, |D). Fig. [2ja,b) 
clearly demonstrate that the spin thermalization is much 
faster than the exciton population decay as was assumed 
in the model. At the lowest temperatures, the population 
decay time is comparable to the laser repetition period, 
which must be taken into account when extracting the 
slow amplitude, A 8 , from these curves. The fast ampli- 
tude, Af , is extracted by comparing the initial part of the 
decay curves in Fig. [21(a) with a local single-exponential 
fit in the time range 2-5 ills (marked by straight lines for 
two of the curves in Fig. [21(a))- 

The extracted values of A s and Af/A Q are plotted 
in Fig. [3](a,b) and compared to the model of Eq. (j4j) 
(red curve). We find a reasonable agreement when 
R = 0.013 ± 0.003 and A = 15.0 ± 1.5 meV. The data in 
Fig. [3](a) is a consequence of thermal equilibrium and es- 
sentially confirms previously reported results^^. The 
observation that the data in Fig. [3](b) follows our model 
is strongly supporting the assumption that the initial de- 
cay (Fig. |2](a)) is caused by spin thermalization between 



the bright and dark exciton states. 

Since the decay curves shown in Fig. [2Ka,b) are in gen- 
eral not single exponential, we make a multi-exponential 
fit, f(t) = ctj exp(— 7j£), to the curves and define the 
characteristic decay rate, 7* = Y^j a j/ J2j Tf- This decay 
rate corresponds to a single-exponential decay preserving 
the initial amplitude and the area under the decay curve. 
The data in Fig. [2Kb) has been fitted using three terms, 
and the resulting characteristic decay time, I/7*, is plot- 
ted in Fig. [31(c). The data ranges from 45 lis at 294 K to 
1.1 ms at 16 K. The fact that the relative change in 7*, 
which depends on Tb and Td, is comparable to the rel- 
ative change in A s , which depends on rB,rad and To^ad, 
indicates that the quantum efficiency of the exciton light 
emission is relatively high. 

The initial thermalization part of Fig. [21(a) is fitted 
using one or two exponential terms, two examples are 
shown in Fig.[2](c,d) for 16 K and 294 K, respectively. As 
exemplified, the lowest temperatures require two terms 
in the fit, while one term is sufficient for the highest tem- 
peratures. The characteristic thermalization time, 1/7^, 
is plotted in Fig. [3](d), and we see that for low temper- 
atures the thermalization time is approximately 100 ns. 
For higher temperatures, the initial random population 
is much closer to thermal equilibrium, which seems to 
slow down the thermalization rate. This slowing down is 
also indicated by the requirement of a bi-exponential fit 
in Fig. [21(c), although we cannot exclude inhomogeneous 
broadening effects on the spin-flip time. In the literature, 
a commonly applied mode l 10 ! 11 for the spin- flip rate is: 
Ted = Fo(N + 1), Tdb = To A 7 ", where To is a charac- 
teristic zero-temperature rate and N = (e A / kT — 1) _1 is 
the number of phonons available at the transition energy, 
A. This model is certainly invalid in our case. Although 
the two-state model in Fig. [U(b) assumes the (unknown) 
splitting 18 between heavy and light holes to be larger 
than /cT, the analysis shows that the model captures the 
main characteristics of the spin- flip dynamics. 

The results can be compared to bright/dark-state spin- 
flip times in other NC systems at low temperatures. For 
InAsi^ and InGaA a 1Q i 13 NCs similar time scales of the 
order of 100 ns have been reported. However, this is not 
unique for all NCs. In CdSeii NCs the spin-flip time is of 
the order of 10 ns, whereas a much faster upper bound of 
200 ps was reported 14 for InP NCs (in this case a constant 
low-temperature spin-flip time was also reported). 

In conclusion, we have measured the spin-flip thermal- 
ization time between bright and dark exciton states in 
silicon nanocrystals and found a constant thermalization 
time of 100 ns below 100 K and counter-intuitively in- 
creasing with temperature above 100 K. This work was 
supported by The Danish Council for Independent Re- 
search I Natural Sciences (FNU) and Technology and 
Production Sciences (FTP, SERBINA project). We are 
grateful to Brian Bech Nielsen for supplying the cryo 
cooler and for useful discussions. 
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